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Abstract: Collisional effects are included in the path-integral formulation 
that was proposed in one of our previous paper for the collisionless Boltz- 
mann gas. In calculating the number of molecules entering a six-dimensional 
phase volume element due to collisions, both the colliding molecules and the 
scattered molecules are allowed to have distributions; thus the calculation 
is done smoothly and no singularities arise. 
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1 Introduction 



In our previous works, we proposed a path-integral approach to the col- 
lisionless Boltzmann gas^Q. It is assumed in the approach that there 
are continuous and discontinuous distribution functions in realistic Boltz- 
mann gases: continuous distribution functions are produced by continuous 
distribution functions that exist previously and discontinuous distribution 
functions are caused by boundary effects. (Boundaries can block and reflect 
molecules in such a way that distribution functions become discontinuous in 
the spatial space as well as in the velocity space.) To treat these two kinds 
of distribution functions at the same time, a different type of distribution 
function, called the solid-angle-average distribution function, is introduced 
as 

f(t,v,v,An) = -^j f(t,r,v)dn, (1) 

where AO represents one of the solid angle ranges in the velocity space 
defined by the investigator and f(t, r, v) is the "ordinary" distribution func- 
tion. By letting each of Af2 be adequately small, the newly employed dis- 
tribution function is capable of describing gas dynamics with any desired 
accuracy. Provided that collisions in a Boltzmann gas can be neglected, the 
solid-angle-average distribution function is found to be 



tf. AO x 1 / / Oo,r ,v,0 ) \cos a\dS 

f{t, r, v, AO) = — / : U ror 

MlJASi |r-r r 

+ AnJ A s 2 Ir-rolV r ° r ' 



(2) 



where, referring to Fig. 1, A5i is an arbitrarily chosen virtual surface within 
the effective cone defined by —Ail at the point r, AS2 stands for all bound- 
ary surfaces within the effective cone, r\ is the local emission rate of boundary 
surface (acting like a surface-like molecular source), ro represents the posi- 
tion of dSo, Qo is the solid angle of the velocity but takes the direction 
of (r — ro), to is the local time defined by to = t ~ \ r ~ r o\/ v i a is the 
angle between the normal of dSo and the vector r — ro, f ct is the contin- 
uous part of the distribution function existing previously, and U ror is the 
path-clearness step function, which is equal to 1 if the path ror is free from 
blocking otherwise it is equal to 0. 

The objective of this paper is to include collisional effects in the path- 
integral formalism. 
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According to the conventional wisdom collisions can be analyzed by the 
method developed by Boltzmann long ago, in which it is understood that 
there is a symmetry between the ways molecules enter and leave a phase vol- 
ume element. Peculiarly enough, this well-accepted understanding includes 
actually hidden fallacies |3|], which can briefly be summarized as follows. In 
terms of studying collisions in a Boltzmann gas, there are two issues that are 
supposedly important. The first one is related to how many collisions take 
place within a phase volume element and during a certain time; the second 
one is related to how the scattered molecules will, after collisions, spread out 
over the velocity space and over the spatial space. These two issues involve 
different physics and have to be formulated differently. If the molecules leav- 
ing a phase volume element is of interest, one needs to take care of only the 
first issue; whereas if the molecules entering a volume element is of inter- 
est, one needs to concern oneself with both the issues aforementioned. This 
imparity simply suggests that the time-reversal symmetry, though indeed 
exists for a single collision between two molecules, cannot play a decisive 
role in studying collective effects of collision. 

In the present paper we formulate the collisional effects partly in an 
unconventional way. In deriving how many molecules make collisions, the 
standard method is employed without much discussion; but, in formulat- 
ing how scattered molecules enter a six-dimensional phase volume element, 
which is an absolute must for the purpose of this paper, a rather different 
and slightly sophisticated approach is introduced. 

In sec. 2, general considerations concerning basic collisional process are 
given. It is pointed out that only the scattering cross section in the center-of- 
mass frame is well defined and can be employed in our studies. Sec. 3 gives 
a formula that describes how a molecule, when moving along its path, will 
survive from collisions. Sec. 4 investigates how molecular collisions create 
molecules that enter a specific phase volume element. In the investigation, 
both the colliding molecules and the scattered molecules are allowed to have 
distributions. (Otherwise, singularities will arise, as Ref. 3 reveals.) Sec. 5 
includes all the collisional effects in a complete path-integral formulation. In 
Sec. 6, approximation methods are introduced to make the new formulation 
more calculable and an application of the method is demonstrated. Sec. 7 
offers a brief summary. 

Throughout this paper, to make our discussion as simple as possible, it 
is assumed that molecules of interest are all identical, but distinguishable, 
perfectly rigid spheres and they move freely when not making collisions. 
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2 General considerations of collision 



Firstly, we recall general features of binary collisions in terms of classical 
mechanics. Consider two molecules: one is called molecule 1 and the other 
molecule 2. Let Vi and v 2 label their respective velocities before the collision. 
The center-of-mass velocity and the velocity of molecule 1 relative to the 
center-of-mass are before the collision 

c = ^( v i+ v 2) and u = ~(vi-v 2 ). (3) 

Similarly, the center-of-mass velocity and the velocity of molecule 1 relative 
to the center-of-mass are after the collision 

c ' = \ W\ + v 2) and u = -(vi - v 2 ). (4) 

The conservation laws of classical mechanics tell us that 

c = c' and |u| = ju'| = u. (5) 

Fig. 2 schematically illustrates the geometrical relationship of these ve- 
locities. Note that, for the collision defined as above the final velocities 
of the two molecules, such as u', v[ and v 2 , cannot be completely deter- 
mined unless the impact parameter of the collision is specified at the very 
beginning]!]]. 

At this point, mention must be made of one misconcept in that the 
usual derivation of the Boltzmann equation gets involved ||]. In an attempt 
to invoke the time-reversal symmetry of molecular collision, the standard 
treatment in textbooks || defines the scattering cross section in the labora- 
tory frame in such a way that 

ct(vi,v 2 v' l5 v 2 )tividv 2 (6) 

represents the number of molecules per unit time (per unit flux of type 1 
molecules incident upon a type 2 molecule) emerging after scattering with 
respective final velocities between v[ and v[ + dv[ and between v' 2 and 
v 2 + dv' 2 . 

If a close look at expression (|6|) is taken, we may find that the value 
of a in it is ill-defined. As Fig. 2b clearly shows, the molecules of type 1, 
namely the ones with the velocities vi after the collisions, will spread out 
over a two-dimensional surface in the velocity space (forming a spherical shell 
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with diameter 2u) rather than over a three-dimensional velocity volume as 
suggested by the definition. Because of this seemingly small fault, the value 
of a actually depends on the size and shape of dv' t and cannot be treated 
as a uniquely defined quantity theoretically and experimentally 

Another type of scattering cross section, which is elaborated nicely in 
textbooks of classical mechanics and suffers from no difficulty, is in terms of 
the relative velocities u and u', as shown in Fig. 3. It is defined in such a 
way that the area element 

dS = a(n u ,)dn u , (7) 

represents the number of molecules per unit time (per unit flux of type 
1 molecules with the relative velocity u incident upon a type 2 molecule) 
emerging after scattering with the final relative velocity u' pointing in a 
direction within the solid angle range dQ u i. Note that the definition (0), 
in which the center-of-mass velocity c or c' becomes irrelevant, makes good 
sense in the center-of-mass frame rather than in the laboratory frame. 

Before finishing this section we turn to discussing how collisions can 
generally affect the solid-angle-average distribution function f(t, r, v, Af2) 
defined by (JlJ). In view of that gas dynamics of the Boltzmann gas develops 
along molecular paths, as shown by (|2|), we believe that collisional effects 
should also be investigated and formulated in terms of molecular paths. 
Fig. 4 illustrates that there are two types of processes. On one hand, 
a molecule that would reach r with the velocity v at time t may suffer 
from a collision and become irrelevant to the distribution function; on the 
other hand, an "irrelevant" molecule may collide with another molecule 
and then become relevant. It should be stressed again that there is no 
symmetry between the two types of processes. For the first type of process, 
we only need to investigate what happens to a single molecule. If a collision 
takes place with it, we know that the molecule will depart from its original 
path, which is sufficient as far as our formulation is concerned. For the 
second type of process, we need to know: (i) how many collisions take 
place within the effective cone; (ii) how the scattered molecules spread out 
over the phase space. As stressed in the introduction, the second issue is 
particularly essential because of that the distribution function is nothing 
but the molecular density per unit phase volume, in other words we must 
concern ourselves with the scattered molecules "around" the phase point 
(r, v), rather than the scattered molecules "at" the phase point (r,v). 

In the next two sections, we will formulate the two processes respectively. 
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Since the collisions are assumed to take place in terms of classical mechanics 
all the calculations can be done without analytical difficulty. 



3 The surviving probability 

Consider a molecule moving along a spatial path where many other molecules 
make their own motions. If P(j) denotes the probability that the molecule 
survives a time r without suffering a collision and w{r)dT denotes the prob- 
ability that the molecule makes a collision between time r and time r + dr, 
we must have a simple relation 

P(t + dr) - P(r) = -P(T)w(T)dT, (8) 

which yields 

1 dP 

Therefore, the surviving probability associated with a molecule moving from 
ro to r with the velocity v can be expressed formally by 

P(r , r; v) = exp(- Jw(r)dr), (10) 

where I represents the path along that the molecule will move if no collision 
takes place. For the Boltzmann gas under consideration, whose molecules 
are assumed to be free from forces except in collisions, the path of a molecule 
is nothing but the segment of straight line linking up the two points. 

The surviving probability defined by (10) can be evaluated by the stan- 



dard approach^. To make this paper complete, we include the evaluation in 
what follows. Suppose that the molecule encounters a molecular beam with 
the velocity vi at the path element dl. In terms of the molecular beam, the 
molecule has the speed 2u = 2|u|, in which u = (v — vi)/2, and it occupies 
the volume with respect to the beam 

2ua(n u ,)dn vil dT, (11) 

where a(Q u i) and Q u i are defined in ([?]) and illustrated by Fig. 3. The 
molecular density of the colliding beam is 

/(r.rj.vOdvi, (12) 
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where r; is the position of the path element dl. Thus, the total collision 
probability can be written as 



wdr = dr j j 2uf(r, ti, vi)cr(Q u ')df2 u 'dvi. 
In terms of (|Hi|), the surviving probability (|To| ) becomes 



(13) 



P(r ,r;v) = exp 



l Jvi JQ„> 



2w(f2 u /)/(r, ri,vi)dQ u i dvidr. 



(14) 



where dr is the time period during that the molecule passes the path element 
dl. 

The formula ( [HP describes how collisions make the number of molecules 
along a certain path decrease. The method employed has nothing partic- 
ularly new in comparison with that employed by the textbook treatment. 
Before changing our subject, one thing worth mentioning. In deriving fll~4]) , 
we had luck not to be concerned with how the scattered molecules spread 
out over the phase space. It is readily understandable that the same luck 
will not be there in the next section. 



4 The creation probability 

We now study the process in which collisions make molecules give contribu- 
tions to the solid-angle-average distribution function f(t, r, v, AO). 

It should be mentioned that in this section, unlike in the last sections, 
v' and v'i represent the velocities of colliding molecules while v and vi 
represent the velocities of scattered molecules. 

The see how the scattered molecules spread out over the six-dimensional 
phase space, we consider a relatively small six-dimensional volume element 
as 

Ar • Av = Ar • v 2 AvAft. (15) 



In (15) Ar is chosen to enclose the point r in /(i, r, v, Ail), Av to enclose the 
speed v in f(t, r, v, Ail); and Aft is just the finite velocity solid-angle-range 
Ail in f(t, r, v, Ail). The discussion below will be focused on molecules that 
really enter, after collisions, this six-dimensional volume element. 

Note that in Fig. 5a each point in the spatial volume Ar defines an 
effective cone, within which physical events may make impact on the distri- 
bution function at the point. This means that the entire effective cone, with 
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respect to the spatial volume element Ar, must be slightly larger than the 
one defined solely by r, as shown by Fig. 5b. Fortunately, we will, in the 
end of formulation, let 

Ar -> (16) 

and thus the actual entire effective cone is only academically larger. On this 
understanding, we will not distinguish between the entire effective cone and 
the effective cone defined solely by the single point r. 

Look at molecular collisions taking place within the effective cone shown 
in Fig. 6. Note that Fig. 6a, while coming to one's mind immediately, 
is not an appropriate picture to manifest the collision process affecting the 
distribution function at the point r since both the colliding molecules and 
scattered molecules in it have no true distributions. (Ref. 3 analyzes the 
situation and brings out that such mental picture will finally lead to singu- 
larities.) In Fig. 6b, the velocity distributions of all colliding molecules and 
scattered molecules are explicitly illustrated. Our task here is to formulate 
the relationship between the distribution functions of colliding molecules 
and the scattered molecules (including their velocity distributions and spa- 
tial distributions). 

We divide the entire effective cone into many individual regions, denoted 
by (Aro)j. It is obvious that within each of the regions collisions can gen- 
erate a certain number of molecules that will finally enter the phase volume 
defined by (|l5|). Let nf denote the number of such molecules. In what 
immediately follows, it is assumed that the generated molecules suffer no 
further collisions. The entire contributions of all collisions to the distribu- 
tion function f(t, r, v, Afl) can then be expressed by 

r^r^A^H^-i^-^nf, (17) 

i 

in which i runs all the divided regions within the effective cone. For later 
use, we wish to rewrite (llTT) as 



rt''''^'^?^ (18) 



The advantage of (18) over (0) is that (Ar)* and (Av) t in (jT|) may be 
chosen to be different for different i as long as the molecular number of nf 
is counted up accordingly. (Of course, all velocity directions of the molecules 
have to be within the solid angle range AO.) 
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Now, consider a small, much smaller than Af2, solid angle range AQq 
at a point ro towards the point r, as shown in Fig. 7a (ro is within the 
effective cone). It is easy to see that if collisions take place at ro, the 
scattered molecules having velocities within A^o will spread out over the 
spatial volume element 

Ar « |r - r \ 2 vAtt At, (19) 

as shown in Fig. 7b. Accordingly, they will spread out over the velocity 
volume element 

Av » v 2 AvAn , (20) 

as shown in Fig. 7c. Since AQq is much smaller than AO, (the latter one 
is finite), a molecule having a velocity within Af^o can be regarded as one 
having a velocity within AO. Thus, by letting (Ar)j in (18) be equal to 



Ar of (|l9j) and letting (Av)i in ( |18|) be equal to Av in (|20f) , expression (|l 
becomes 

nt,r,v,An) = JL / l im _ , pd(r0) f 2 r A ° o , (21) 

where the integral is over the entire effective cone defined by r and —Ail 
and p cl is the local density (per unit spatial volume) of the molecules that 
are "emitted" from ro due to collisions and finally enter the speed range Av 
and the solid-angle range Af^o during the time At. 

To determine the density p d , we have two tasks. One is to derive the 
collision rate at ro and the other is to derive what fraction of the scattered 
molecules emerge with velocities within the range AvAQq. The first task 
can be accomplished in a well-known way while the second one cannot. 

As discussed in the last section, a specific molecule with the initial ve- 
locity v' that collides with a beam of molecules with the initial velocity 
occupies a volume with respect to the beam 

2uAta(Q u )dQ u , (22) 

where Q u is the solid angle of the scattered relative velocity u = (v — vi)/2. 
The number of "such specific" molecules within dv'dro is 

f{t ,r o y)dv'dr . (23) 

The density of the molecules with is characterized by 

/(i ,r ,vi)c/vi. (24) 
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Therefore, the number of all collisions within the spatial volume dro within 
the time At is 



dr J dv' J dv[J dn u /(v')/(vi)2uo-(ft u )Ai. (25) 



We now to evaluate the probability that the molecules expressed by fl2q ) 
enter the velocity range AvAQq. Note that the integration of dv'dv'i is 
carried out in the laboratory frame while the integration of cZO u is in the 
center-of-mass frame, which makes the evaluation quite difficult. For this 
reason, we make the integration conversion as 

J dV f dv[--- = J dc' J dfl u , J u 2 du\\J\\ • • • , (26) 

where u = u' is understood and ||J|| represents the Jacobian between the 
center-of-mass frame and the laboratory frame 

II Jll " d(V ' yi) (27) 

Equation (Q) tells us that the Jacobian is equal to 8. 

By making use of ( |25| ) and (^6|), the distribution function f cl expressed 
by (|2l|) becomes 



f d (t,r,v,An) rs — / dr f dc' f dQ u , f f u 2 dn u du 

Ail J-An J J J JAvAtto , , 

IIJII ( 2 °) 
2ua(tt VL )f(t ,r ,c' - u')f(t , r , c' + u ). 



<■'■'> Ac r - r'| 2 Afi 



In regard to expression (|2^), some observations are made. As mentioned 
in Sec. 2, if two molecular beams with definite velocities v' and v^ col- 
lide with each other the scattered molecules will spread out only over a 
two-dimensional spherical surface in the velocity space, which implies that 
difficulty arises if the velocity distributions of scattered molecules are of con- 
cern. Whereas, in this expression, all the colliding molecules are allowed to 
have distributions, the value of u' = u is allowed to vary and therefore the 
scattered molecules explicitly spread out over the velocity space (as well as 
over the spatial space). Furthermore, by using the notation 



(29) 

AvAUo 
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we have ensured that only the scattered molecules of relevance are taken 
into account. 

In Fig. 8, which is drawn for scattered molecules in the velocity space, 
we are concerned only with molecules that finally enter the range AvAQq. 
Allowing u to vary a little bit, we may let the scattered molecules fill out 
the velocity range. Namely, we have 

u 2 dn u du(---) ^v 2 AvAn (---), (30) 

where (• ■ ■) represents other factors that have been treated as constants in 
terms of the infinitesimally small range of AvAQq. 

Inserting (|30|) into PH ) and taking the limits AQq — ► and Av — > 0, we 
finally arrive at 

f d (t,r,v,Afl) = -L. [ dr f dc' f dn u , 

VAll J-AQ J J ^g-g 

— 2ua(Q u )f{t , r , c' - u')f(t , r , c' + u'), 
|r — ro| 

where to = t — \r — tq\/v and the integration 

dQ u , ■ ■ ■ (32) 

is over the entire solid angle (0 — > Air). Note that u, u' and u in the 
integrand of ( |3l|) have to be determined skillfully. First use c = c' and 
v = v(r — ro)/|r — ro| to determine u, then use u = |u| and Q u i to determine 
u', as shown in Fig. 9. 

We have directly formulated the contribution to the solid-angle-average 
distribution function from collisions. It should be noted that the formula- 
tion can be done only under the condition that the velocity solid-angle range 
AO is finite: if both AO and AO,q in the formulation were assumed to be 
infinitesimally small, the limiting processes concerning the two quantities 
would not be in harmony with each other. This shows again that the intro- 
duction of the solid-angle-average distribution function is a must to the gas 
dynamics of Boltzmann gas. 



5 Complete formulation 

The complete formulation for the collisional Boltzmann gas is now in order. 
In Fig. 10, we have depicted a piece of boundary and some collisions taking 
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place within the effective cone. As said before, all these events can directly 
contribute to the solid-angle-average distribution function. 

We then use the following sum to represent the total distribution function 

f(t, r, v, Aft) = + + (33) 

where f^, fuj) and f(ia) represent the contributions from the existing con- 
tinuous distribution function, from the piece of boundary and from the col- 
lisions respectively. For simplicity, no other types of distribution functions 
are assumed to exist within the effective cone. 

As has been illustrated in Sec. 3, a molecule that makes its motion 
towards the point r may suffer a collision with other molecules. The involved 



surviving probability P has been defined by (14). By taking the probability 
into account, the first term in (||) becomes 

t u Ar>\ 1 f f ct (to,r ,v,n )\cosa\dS „, > foA . 

/ ( „(t, r,t,, AO) = ^ P(r ,r,v ), (34) 

where vo in P take the value of v and points to the direction of (r — ro). 

By taking the same effect into account, the second term in (|2|) can be 
expressed by 

f(ii)(t,r,v, An) = — / : |2-g P(r 0) r,v ). 35 

AS I Jas 2 \ r — r o| v 

A rather detailed discussion on the molecular emission rate r\ has been 
included in Ref. 1. Here, we content ourselves with pointing out that the 
rate r\ satisfies the normalization condition at the surface element dSo. If no 
molecular absorption and production by the surface element are assumed, 
the following expression holds 

J rj(t, r, vi, 0,i)dUidvi = J V2f(t,r,V2,Q.2)\cos0\v2d{l2dv2, (36) 

in which 9 is the angle between the velocity V2 and the normal of dSo, fl± 
points to an outward direction of dSo and SI2 points to an inward direction of 
dSo. The concrete relation between rj(t, r, v) and f(t,r, v) must, of course, 
be ultimately determined by experimental data((|. 

In obtaining (^Tj) for the distribution function created by collisions, fur- 
ther collisions were excluded. To include possible further collisions, the 
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contribution expressed by ( J3l| ) needs to be modified as 
f(in){t,r,v, AU) = J An dr ° / dc ' J d ^ u ' 



\J\ 



|r - r | 2 



2w(Q u /)/(i ,r ,c' - u')/(io,ro,c' + u')P(r , r, v ) 



(37) 



where the integration of dro is over the entire effective cone defined by r 
and —Ail, |vo| = v and takes the direction of (r — ro), u' is defined by 
u = |c — voj and Q u /. 

In these formulas, the probability P should be set to be zero at the very 
beginning if there is physical blocking along the path Yqt. 

Equations (|33|)-(37) constitute a complete set of integral equations that 



describe the collisional Boltzmann gas defined in this paper. The formulation 
proves in a theoretical way an obvious intuition that the distribution function 
at a spatial point can directly be affected by physical events taking place 
at other, even remote, points in view of the fact that a molecule can freely 
pass any distance in a certain probability. In this sense, the picture here is 
more "kinetic" than that associated with the Boltzmann equation, in which 
physical events have to make their influence region by region (like what 
happen in a continuous medium). 

Another comment is about the famous H-theorem. If the involved 
distribution function is initially nonuniform in the spatial space and non- 
Maxwellian in the velocity space, the resultant distribution function given 
by the formalism will approach the uniform Maxwellian. Though such ex- 
plicit proof has not been achieved yet, we believe that this must be the case 
by noticing a general discussion stating that as long as a statistical process 
is a Markoffian one the .ff -theorem must hold trueQ. 



6 Approximation and application 

Although the formulation offered in the previous section is formally com- 
plete, there still exist difficulties that hinder one from performing calculation 
for a realistic gas. Unlike the solution for the collisionless Boltzmann gas, 
given by (||), the equation system in the last section, namely (|33|)-(j3"7|), is an 
integral-equation set. Without knowing the entire history of the distribution 
function f(t), the integrals in the system cannot be evaluated accurately. 

Fortunately, there are situations for which adequate approximations can 
be introduced and meaningful results can be obtained. In what follows, we 
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first deal with weakly collisional gases and then give some discussion on how 
the consideration can apply to more general cases. 

If the density of a Boltzmann gas is relatively low, by which we imply that 
the mean free path of molecules is not too short comparing with the length 
scale of the system or that the mean free time is not too short comparing 
with the time scale of the phenomena of interest, we may apply the following 
iterating procedure to calculate the distribution function. 

Firstly, we assume that the system can be treated as a collisionless Boltz- 
mann gas and the collisionless formulation can directly applied. Namely, we 
have the zeroth-order solution 



AO] = J_ 

AQ j as 



o 



f (to,r ,v,n )\ cos a\dS ( 
|r - r | 2 
rt(t ,r Q ,v,Sl )dS 



(«) 



j 3 AQ 



AS 2 



r o| 



(38) 
(39) 



Then, we can construct the first-order distribution function by inserting 
the zeroth-order solution into all right sides of the equations (p^)-(]3^), which 
yields 



[i] 



1 



(*) AO 



f ct (t ,r ,v,Q ) \ cos a\dS ( 



1 

v 3 AO Jas 2 



|r - r | z 
7](t ,r 0l v,tt )dSo 
|r - r | 2 



P [OJ (ro,r,v ) 



P[°](r ,r,v ). 



(40) 
(41) 



and 



J (Hi) 



1 



drQ I dc' / dQ u i- 
J J 



\J\ 



vAQ J -An J J |r — r | 
2w(Q u )/[°] (t , r , c' - u')/ [0] (t 0j r 0j c' + u')P [0] (r Q , r, v ) 



(42) 



In all the first-order formulas, the surviving probability is defined as 



p[0]( 



ro,r; v J 



exp 



Iff 2ua(n u ,)fW(T 1 r l ,v 1 )dn u ,dv 1 dT. 



(43) 



In equations (E^) and (03|), is the total zeroth-order distribution function, 



namely /M = /g + fg. 

Along this line, we can formulate higher-order solutions for dilute gases. 
If the gas of interest is rather dense, the approximation method presented 
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above may not work effectively. One wishes, however, to point out that for 
the regions near boundaries, where the distribution function suffers from 
most irregularities and collisions between molecules have no enough time 
to erase such irregularities, the introduced method should still make sense. 
In view of this, it is expected that a hybrid method will be developed, in 
which the approach here and other effective approaches, such as the ordinary 
fluid theory, can be combined into one scheme so that more practical gases 
become treatable. 

To illustrate the application of our approximation scheme, we investigate 
a gas leaking out of box through a small hole. For simplicity, we assume, 
referring Fig. 11, that the zeroth-order solution of the leaking gas is con- 
fined to a "one-dimensional thin pipe" (shaded in the figure), which can be 
expressed by 

f[0] = f Co exp [-mvl/(2KT)] (inside the pipe) ^ 
| (outside the pipe) ^ ' 

and then we try to determine the collisional effects of the distribution func- 
tion. 

Note that the distribution function expressed by (EJj) is kind of special so 
that the formula (42) should be slightly modified. For this purpose, we write 
the differential collision probability as (the subindex x of v x is suppressed) 

[f(v')ASdx dv'} [f(v[)dv[] [2ua(n u )dn u ], (45) 

where AS is the cross area of the pipe. By making the variable transforma- 
tion, we obtain 

dv' / dv[---= dc' / du'\\J\\--- 



oo 

2 m.a/„.2\ 



dc' / u 2 duA(u 2 



In a way similar to that has been presented in the last section, we finally 
arrive at 

oo 

^ A ^] = ^ / dxo [dc'^^- 2 fM(c' + u')fM(c'-u'), m 
v{A0)i J J u r - r r 

-(Afl), o 

where (A#)j is the polar angle range set by the investigator (the azimuthal 
angle range is irrelevant in the case). 
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The formula ([!(]) can be calculated easily with a computer. Referring to 
Fig. lib, we set 

i 1 m i 

v = l rj_ = 1, r-^= = 1, 



2kT 



let (A0)j be the interval 



0.4vr - 0.52 + — , 0.4vr - 0.5 + — 
50 50 



and notice <r(f2 u ) is constant [0]. The numerical results are listed as the 
following: (normalized by /[(A0)o]) 



f[(A6) ] = 1.00000e + 00 
/[(A6») 5 ] = 5.23910e - 01 
/[(A0)i O ] = 2.01786e - 01 
/[(A6»)i 5 ] = 5.03200e - 02 
/[(A0) 2O ] = 4.76340e - 03 
/[(A0) 25 ] = 5.11013e - 05. 



(47) 



7 Summary 

In this paper, we have proposed a complete mathematical scheme to deal 
with the Boltzmann gas. The scheme has many new features. In addition 
to those given in Ref. 1, some related to treating collisional effects are the 
following. 

Firstly, collisional effects are investigated in the full velocity-and-position 
space. In particular, a six-dimensional volume element is explicitly defined 
and a calculation concerning molecules entering the volume element is di- 
rectly performed. 

Secondly, both the colliding molecules and scattered molecules are al- 
lowed to have distributions. In other words, we consider the full and collec- 
tive behavior of collisions, in which the time-reversal symmetry existing for 
a collision of two molecules plays almost no role. 

Thirdly, the treatment in this approach is consistent with the previous 
approach to the collisionless Boltzmann gas in the sense that all the formulas 
are given in terms of what happen along molecular paths. 

Finally, the resultant formulas of this approach are, in many practical 
situations, calculable by means of today's computer. 

It is believed that this approach will be developed further so that a better 
understanding of complicated fluid phenomena can be achieved. 
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Figure captions 

1. A physical surface and a virtual surface within the effective cone de- 
fined by r and Afi. 

2. A collision between two molecules, (a) The molecular velocities before 
the collision, (b) The molecular velocities after the collision. 

3. The scattering cross section in the center-of-mass frame, (a) Solid an- 
gles and relative velocities, (b) The relation between the cross section 
and solid angle range. 

4. Two types of collision processes. 

5. Effective cones, (a) For a single spatial point, (b) For a given spatial 
volume. 

6. (a) A mental picture in which two molecular beams collide with each 
other, (b) A mental picture in which both colliding molecules and 
scattered molecules have distributions. 

7. (a) A solid angle range AO,q towards the point r. (b) The distribution 
of scattered molecules in the spatial space, (c) The distribution of 
scattered molecules in the velocity space. 

8. The velocity distribution of scattered molecules in the center-of-mass 
frame and in the laboratory frame. 

9. Relations between various essential vectors in the formulation. 

10. Contribution to the solid-angle-average distribution function from dif- 
ferent sources. 

11. A gas leaking out of a container through a small hole. 
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Figure 1 




Figure 2 




Figure 3 
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Figure 4 
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Figure 7 
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Figure 10 
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